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Abstract 

A  unified,  three-dimensional,  steady-state  numerical  mass-transfer  single  cell  model  for  polymer  electrolyte  membrane  fuel  cell  (PEMFC) 
was  developed.  The  modeled  fuel  cell  uses  metal  foam  in  the  flow-field  of  the  bipolar/end  plates  instead  of  the  conventionally  used  rectangular 
channels.  Transport  equations  formulated  under  the  PEMFC  conditions  were  solved  using  the  commercial  computational  fluid  dynamics 
software  Fluent®  6.0  with  Gambit®  2.0  as  pre-processor.  Simulations  were  performed  for  different  permeability  levels  of  the  metal  foam  in 
the  flow-field.  Results  showed  a  significant  effect  of  permeability  of  the  metal  foam  on  the  performance  of  the  fuel  cell.  For  example:  at 
10-6  m2  permeability  of  metal  foam  the  value  of  average  current  density  was  5943  A/m2  while  at  10~u  m2  permeability,  the  average  current 
density  was  8325  A/m2.  The  average  current  density  value  for  the  multi-parallel  flow-field  channel  design  (channel  width  =  0.0625  in., 
channel  depth  =  0.0625  in.  and  land  width  =  0.0625  in.),  which  corresponded  to  an  equivalent  permeability  value  of  4.4  x  10-8  m2  was 
7019  A/m2.  This  value  for  the  porous  configuration  with  same  permeability  and  under  similar  conditions  of  temperature,  pressure  and 
reactants  flow  rate  was  slightly  lower  at  6794  A/m2.  The  trend  indicated  that  decreasing  the  permeability  of  the  flow-field  results  in  better 
performance  from  the  cell.  However,  the  permeability  of  the  channel  design  can  not  be  decreased  below  the  value  of  around  10-8  m2,  due  to 
difficulty  in  machining  thinner  channels.  Consequently,  the  use  of  metal  foam  flow-field  is  proposed  in  the  bipolar/end  plate.  The  developed 
model  offers  fuel  cell  developers  a  scope  for  improvement  of  the  bipolar/end  plates  in  the  fuel  cell,  by  switching  over  to  the  metal  foam  flow- 
field  concept. 
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1.  Introduction 

Fuel  cell  vehicles  (FCVs)  may  be  the  only  option  to  meet 
new  EPA  emission  standards,  according  to  the  National 
Research  Council  “Review  of  the  Research  Program  of 
the  Partnership  for  a  New  Generation  of  Vehicles  (PNGV) 
Sixth  Report.”  Despite  significant  technological  progress 
made  with  PNGV  prototypes  unveiled  this  year  (all  hybrid 
AFVs),  the  report  states  that  production  models  slated  for 
debut  in  2004  will  not  meet  new  EPA  emission  standards  [1]. 
Fuel  cell  vehicles  have  the  greatest  potential  to  meet  these 
standards  and  fuel-efficiency  requirements.  With  the  launch 
of  FreedomCAR  (Cooperative  Automotive  Research)  pro¬ 
gram  by  US  administration  in  early  2002,  the  focus  is  now  on 
mass  production  and  commercialization  of  hydrogen  pow¬ 
ered  FCVs.  The  FreedomCAR  program  will  focus  on  devel¬ 
opment  of  fuel  cell  technologies  and  efficient  production  of 
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hydrogen  from  domestic  renewable  sources.  The  program 
replaces  the  PNGV  initiative  to  make  fuel  efficient  vehicles 
launched  by  the  US  administration  in  1993. 

It  is  clear  that  the  fuel  cell  technology  has  the  right 
potential  to  replace  the  conventionally  used  internal  com¬ 
bustion  (IC)  engine  technology  in  the  automobiles.  Among 
the  different  types  of  fuel  cells  that  exist,  the  polymer 
electrolyte  membrane  fuel  cell  (PEMFC)  is  seen  as  a  system 
of  choice  for  automobile  applications  due  to  its  environ¬ 
mental  friendly  nature  and  ability  to  deliver  high  power 
density  while  operating  at  lower  temperatures  [2] .  With  so 
much  of  progress  in  fuel  cell  technology  in  last  decade,  one 
of  the  main  hindrances  in  the  onset  of  fuel  cell  commercia¬ 
lization  program  is  the  high  cost  of  the  fuel  cell  system.  The 
cost  of  fuel  cell  stack  is  currently  at  US$  500/kW  while  the 
2004  PNGV  goal  is  US$  50/kW  [3].  Therefore,  the  main 
focus  for  the  fuel  cell  developers  is  to  improve  the  perfor¬ 
mance  of  different  components  of  fuel  cell  stack  by  choosing 
high  performance  materials  and  thereby  increasing  the  cell 
efficiency. 
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One  such  component  is  the  fuel  cell  bipolar/end  plates. 
These  are  one  of  the  most  important  and  costliest  compo¬ 
nents  of  the  PEMFC  stack  and  accounts  for  more  than  80% 
of  the  total  weight  of  the  stack  [4].  These  plates  act  as 
separator  plates  in  between  the  different  cells  of  the  stack, 
thereby  keeping  the  oxidant  and  fuel  from  coming  in  direct 
contact.  They  support  intricate  gas  flow-field  patterns  that 
help  to  direct  fuel  and  oxidant  stream  to  individual  cells  in 
the  stack.  They  also  act  as  current  collectors  and  provide  a 
series  of  electrical  connections  in  between  individual  cells. 
The  overall  efficiency  of  PEMFC  stack  depends  on  the 
performance  of  these  plates  in  the  fuel  cell  environment. 
Although  several  experiments  and  model  have  been  devel¬ 
oped  for  the  electrochemical  process  in  the  fuel  cell,  these 
bipolar/end  plates  have  received  little  attention.  The  present 
work  focuses  on  the  improvement  in  the  performance  of  the 
fuel  cell  by  optimizing  the  gas  flow-field  in  the  bipolar/end 
plate. 

The  current  design  for  the  flow-field  uses  the  machined 
rectangular  flow  channels.  Since  the  flow-field  helps  to 
distribute  the  reactant  gases  on  the  surface  of  the  electrode, 
and  remove  the  by-product  water,  the  dimensions  and 
shape  of  the  channels  in  the  flow-field  will  affect  the 
performance  of  the  fuel  cell.  Studies  have  been  done  for 
different  channel  dimensions,  and  it  was  seen  that  smaller 
the  channel  dimensions  (channel  width,  channel  depth  and 
land  width)  and  more  the  number  of  channels  in  the  flow- 
field,  better  is  the  fuel  cell  performance  [5,6].  For  the 
channel  type  of  flow-field,  Kumar  and  Reddy  showed  that 
optimum  dimensions  exists  for  the  channels  that  would 
yield  maximum  performance  from  the  fuel  cell  [7].  How¬ 
ever,  in  general,  lower  is  the  permeability  of  the  flow- field 
better  would  be  its  performance.  The  low  permeability 
would  result  in  increased  pressure  drop  across  the  flow- 
field.  This  will  make  the  transfer  of  reactant  gases  towards 
the  electrode-membrane  reaction  interface  from  diffusion 
to  forced  convection  type.  However,  the  permeability  value 
of  the  flow-field  can  not  be  decreased  beyond  10-8  m2  in 
case  of  channel  design  due  to  difficulty  in  machining  the 
thin  cross-section  channels. 

So  there  is  a  need  to  look  for  alternative  technologies 
which  can  produce  thin  cross-section  channels,  or  at  least 
given  performance  equivalent  to  that  of  thin  channels.  The 
use  of  foam  materials  with  an  open-porous  structure  in  this 
regard  seems  promising.  Metal  foams  are  a  new  kind  of 
material  which  as  on  now  have  not  been  characterized 
properly,  but  do  have  some  alluring  properties.  They  are 
light  and  stiff,  have  good  energy  absorbing  characteristics, 
and  good  heat  transfer  properties.  This  work  therefore 
focuses  on  the  development  of  a  three-dimensional 
steady- state  numerical  mass-transfer  model  for  PEMFC 
with  metal  foam  in  the  flow-field  of  the  bipolar/end  plates. 
The  model  would  then  be  extended  for  the  conventional 
channel  design  flow-field.  The  results  obtained  from  the 
simulations  in  both  the  cases  would  be  studied  and  com¬ 
pared. 


2.  Model  development 

In  this  research  work,  a  three-dimensional  steady-state 
numerical  mass-transfer  unified  model  for  PEMFC  was 
developed  to  predict  the  cell  performance.  Numerical  pre¬ 
dictions  of  current  density  as  a  function  of  permeability  levels 
of  the  flow-field  would  be  made.  This  would  help  in  finding 
the  best  permeability  values  for  the  metal  foams  which  need 
to  be  used  in  the  flow-field.  The  detailed  step-wise  develop¬ 
ment  of  the  model  is  discussed  in  following  sections. 

2.7.  Problem  domain 

The  simulation  domain  consists  of  two  flow-fields  with 
metal  foams  corresponding  to  cathode  and  anode  sides 
separated  by  membrane  electrode  assembly  (ME A).  The 
MEA  consists  of  anode  and  cathode  porous  diffusion  layer 
and  catalyst  layer,  on  each  side  of  membrane.  The  active 
area  of  the  domain  was  a  square  of  side  7.62  cm.  Fig.  1 
shows  a  detailed  view  of  the  model  domain  ( X-Z  section), 
which  is  composed  of  a  cathode  metal  foam  flow-field,  a 
cathode  diffusion  layer,  a  cathode  catalyst  layer,  an  anode 
catalyst  layer,  an  anode  diffusion  layer,  and  an  anode  metal 
foam  flow-field.  Each  of  the  layers  is  marked  as  separate 
zone.  The  diffusion  layers  were  essentially  the  electrodes 
made  with  porous  carbon  cloth.  These  layers  help  to  diffuse 
the  reactant  gases  from  the  bipolar/end  plate  flow  channels 
onto  the  reaction  catalyst  layer  on  each  sides  of  the  mem¬ 
brane.  It  was  assumed  that  these  layers  have  an  isotropic 
permeability  and  the  value  is  10  12  m2  [8].  The  catalyst 
layer  is  where  all  the  electrochemical  reactions  take  place. 
Four  species,  viz.  hydrogen,  oxygen,  nitrogen,  and  water 
vapor  were  considered  in  the  modeling  work.  The  model 
works  by  tracking  the  fluxes  of  all  the  reactants  gases  along 
with  the  by-product  water  which  is  formed  due  to  the 
electrochemical  reaction.  Consider,  for  example,  the  flux 
for  hydrogen  species.  The  hydrogen  from  the  anode  flow- 
field  is  transported  through  the  anode  diffusion  layer  toward 
the  anode  catalyst  layer.  The  catalyst  action  dissociates  the 
hydrogen  molecule  to  protons  and  electrons  according  to  the 
reaction  H2  — >  2H+  +  2e_.  The  protons  generated  above 
will  transfer  through  the  membrane  (protonic  conductor) 
onto  the  cathode  side,  wherein  they  combine  with  oxygen  to 
generate  water.  So  there  is  a  consumption  of  hydrogen  and 
oxygen  species  on  the  anode  and  cathode  side,  respectively, 
and  generation  of  water  species  on  the  cathode  side.  Simi¬ 
larly  other  species  are  tracked  in  the  domain. 

2.2.  Model  assumptions 

The  following  assumptions  were  used  in  developing  the 
unified  PEMFC  model. 

(1)  Steady-state  and  stationary  conditions  exists  in  the 
single  cell  stack.  Also,  the  effect  of  gravity  was 
neglected. 
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Fig.  1.  Schematic  of  polymer  electrolyte  membrane  fuel  cell  cross-section  showing  different  zones  and  species  transport  across  the  zones.  The  net  water  flux 
is  sum  of:  (1)  electro-osmotic  effect,  (2)  diffusion  effect  and  (3)  permeability  effect. 


(2)  Based  on  the  Reynolds  number  calculation  the  flow  in 
the  fuel  cell  domain  is  laminar.  Hence,  all  the  transport 
equations  were  formulated  for  laminar  behavior. 

(3)  All  the  porous  zones  in  the  fuel  cell  domain  were 
assumed  to  be  isotropic. 

(4)  Four  reactant  species,  viz.  H2,  02,  N2  and  H20  (vapor) 
were  considered. 

(5)  The  volume  of  the  liquid  H20  was  assumed  to  be 
negligible  in  the  domain. 

(6)  Since  the  density  of  the  gaseous  mixture  varies  from 
location  to  location  in  the  domain,  compressible  gas 
technique  was  used  to  determine  the  density  of  the  gas 
mixture. 

(7)  Gas  mixture  viscosity  was  calculated  using  mass- 
weighted-mixing-law. 

(8)  The  species  binary  diffusion  coefficient  was  calculated 
using  kinetic  theory  of  gases. 

(9)  The  oxide  film  thickness  on  the  bipolar  plate  was 
assumed  to  be  constant  and  that  the  contact  resistance 
of  the  bipolar/end  plates  with  the  electrodes  was  a 
function  of  compaction  pressure  of  the  stack. 

2.3.  Governing  equations 

The  governing  equations  for  the  model  can  be  divided  into 

three  different  categories,  viz.  (a)  transport  equations,  (b) 

thermodynamic  equations,  and  (c)  membrane  properties 

equations. 


2.3.1.  Transport  equations 

These  include  equations  of  continuity,  momentum  con¬ 
servation  equation,  and  species  transport  equations. 


2. 3. 1.1.  Equation  of  continuity.  The  equations  of  continuity 
or  mass  conservation  for  steady-state  system  can  be  written 
as 


V(pv)  =  Sm  (1) 

where  p  is  the  density  of  the  gaseous  mixture,  v  the  velocity 
vector,  and  Sm  the  source  term  for  the  continuous  phase.  The 
value  of  Sm  is  zero  for  all  other  zones  in  fuel  cell  domain, 
except  that  in  catalyst  layer  for  both  anode  and  cathode 
sides.  The  source  term  in  the  anode  and  cathode  catalyst 
layer  is  given  as 


^m,a  —  ^H2  T~  5\v, a 

^m,c  —  ^02  T~  $w,c 


(2) 

(3) 


where  the  subscripts  ‘a’  and  ‘c’  refer  to  the  anode  and 
cathode  side,  respectively,  and  Sx  (x  =  H2,  02,  w  (water)) 
represent  the  source  terms  for  species  x  in  the  catalyst  layer. 
The  density  (p)  of  the  mixture  is  calculated  using  the 
compressible  gas  technique  [9]  and  is  given  for  multi- 
component  system  as 


P  = 


op 


RTj2,m,/M, 


(4) 
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where  Pop  is  the  operating  pressure,  P  the  local  relative  (or 
gauge)  pressure,  mt  the  mass-fraction  of  species  /,  and  Mt  the 
molecular  weight  of  species  i. 

2. 3. 1.2.  Momentum  conservation  equations.  The  momentum 
conservation  equations  or  the  Navier-Stokes  equations  for 
steady-state  system  can  be  described  as 

V  (pvv  +  P)  +sp  =  5m)kv  (5) 

where  Sp  is  the  momentum  source  term.  For  small  velocities 
and  transport  in  porous  media,  the  momentum  source  term 
Sp  is  given  by  Darcy’s  Law: 

Sp  =  -(j)?  «S) 

where  /i  is  the  viscosity  of  the  gaseous  mixture  given  by 
mass- weighted  mixing  law: 

r  =  (?) 

i 


s 


w,a  — 


F 


I(x,y)Mu2oAC\ 


So2 


i&y) 

4  F 


Mo2ACv 


(15) 

(16) 


S 


w,c 


1  +  2oc(x,y) 
2 F 


7(x,y)MH2o^cv 


(17) 


where  I(x ,  y)  is  the  local  current  density,  F  the  Faradays 
constant,  a(x,  y)  the  local  net  water  transfer  coefficient  per 
proton  and  ACv  the  specific  surface  area  of  control  volume 
element  in  the  domain. 


2.3.2.  Thermodynamic  equations 

The  various  reactions  for  a  PEMFC  fed  with  hydrogen 
containing  anode  gas  and  oxygen  containing  cathode  gas  are 


anodic  :  H2  — >  2H+  +  2e 

(18) 

cathodic  :  ^02  +  2H+  +  2e_  — >  H20 

(19) 

and  is  the  permeability  of  the  medium.  The  value  of  /?  was 
taken  to  be  10“ 12  m2  for  the  electrodes  and  catalyst  layer. 
The  permeability,  ft  for  the  metal  foam  flow-field  in  the 
bipolar/end  plates  was  varied  from  10-6  to  10-12  m2. 


2. 3. 1.3.  Species  transport  equations.  These  are  written  for 
each  of  the  species  H2,  02,  N2,  and  H20  (vapor).  On  the 
anode  side,  the  species  transport  equations  are 


V(pvmH2)  =  -V  •  7h2  +  SU2  (8) 

V(pvww  a)  V  •  7W, a  T~  Sw, a  (9) 

and  for  the  cathode  side  the  species  transport  equations  of 
02,  N2,  and  H20  are 

V(pvm02)  =  -V  •  Jo2  +  S0l  (10) 

V(pvmN2)  =  -  V  -7n2  (11) 

V(pv/tzw  c)  —  V  *  7W  c  H-  Syj^ c  (12) 


where  J,  is  the  diffusion  flux  of  species  i. 

The  diffusion  flux  is  given  by  Maxwell  relationship 
[9]: 

N- 1 

Ji  =  —'y^pDjjVnij  (13) 

7=1 


where  Dtj  is  the  binary  diffusion  coefficient  of  species  i  in 
species  /,  and  N  the  total  number  of  species  in  the  mixture. 
The  binary  diffusion  coefficients  (DUj)  were  calculated  using 
the  kinetic  theory  of  gases  as  described  in  Fluent®  user 
manual  [9].  The  source  terms  for  each  of  the  species 
transport  equations  exist  only  in  their  respective  catalyst 
layer  and  are  given  by 


% 


IF 


^h2^cv 


(14) 


net  reaction  :  H2  +  ^02  — »  H20.  (20) 

The  Nernst  equation  for  the  H2/02  fuel  cell,  using  literature 
values  of  the  standard- state  entropy  can  be  written  as  [10]: 

E  =  1.229  -  (8.5  x  10“4)(r  -  298.15) 

+  (4.308  x  10-5)r(lnQ4)  +  ±  In (p^))  (21) 

where  E  (volts)  is  the  ideal  equilibrium  potential,  T  (K)  the 
cell  temperature  and  P*  (atm)  is  the  partial  pressures  of 
species  i  at  the  reaction  interface.  However,  due  to  the 
irreversible  losses  in  a  practical  fuel  cell  the  actual  cell 
potential  is  decreased  from  its  equilibrium  potential.  These 
losses  include  activation  polarization  (rjact),  ohmic  polariza¬ 
tion  (^ohmX  and  concentration  polarization  (qconc).  The 
actual  cell  voltage  (V)  is  then,  theoretical  cell  voltage  ( E ) 
minus  the  polarization  losses. 

^7  ^  hart  h  ohm  hconc'  (22) 

Activation  polarization  is  present  when  the  rate  of  an 
electrochemical  reaction  at  an  electrode  surface  is  controlled 
by  sluggish  electrode  kinetics.  The  loss  arises  due  to  the 
activation  barrier  that  must  be  overcome  by  the  reacting 
species.  This  is  given  by  the  equation  [2]: 

where  I0  is  the  exchange  current  density  (A/m2). 

The  ohmic  polarization  occurs  because  of  the  resistance  to 
flow  of  electrons  in  cell  components  and  resistance  to  the 
flow  of  ions  in  the  electrolyte,  and  is  given  by  [2]: 

fl0hm=I(x,y)Rceii(x,y)  (24) 

flcell  (*,  y)  =  RCR  (x,  y )  +  Rmemb  (x,  y)  (25) 

where  Rqr(x ,  y)  is  the  contact  resistance  of  the  stack  and 
7?memb  the  resistance  offered  by  the  electrolyte  membrane  to 
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the  flow  of  H+  ions.  The  contact  resistance  depends  upon  the 
type  of  material  for  the  metal  foam  in  the  flow-field  of  the 
bipolar/end  plate  and  compaction  force  of  the  stack.  In  this 
model,  SS-316  metal  foam  was  used.  A  relationship  between 
compaction  force  and  contact  resistance  for  SS-316  was 
developed  by  curve  fitting  the  data  from  Hodgson  et  al.  [1 1]. 
The  equation  for  contact  resistance  for  SS-316  is  given  by: 

R{ §>-316  =  2.62  F“ka9231  (26) 

where  Rq |"316  is  the  contact  resistance  (Qm2)  and  Fstk  the 
compaction  pressure  of  the  stack  in  N/m2.  The  high  contact 
resistance  in  case  of  using  metallic  components  is  due  to  the 
formation  of  thin  oxide  layer  on  the  surface.  Hence,  these 
contact  losses  are  considered  in  our  modeling  calculations. 
The  resistance  to  flow  of  ions  (H+)  in  the  electrolyte 
membrane  depends  on  the  proton  conductivity  of  the  mem¬ 
brane  ((7m)  and  the  thickness  of  the  membrane  (tm),  and  is 
given  by  the  expression: 

flmemb  =  r  V  (27) 


It  may  be  noted  that  the  electrolyte  membrane  is  a 
conductor  of  protons  (H+)  and  is  an  insulator  for  the 
electrons.  Consequently  only  the  protons  will  be  passing 
the  electrolyte  membrane  from  the  anode  to  cathode  side, 
while  the  electrons  will  follow  the  outer  circuit  generating 
current. 

Concentration  polarization  arises  due  to  the  consumption 
of  the  reactant  gases  at  the  reaction  site  and  the  inability  of 
the  surrounding  material  to  maintain  the  initial  concentra¬ 
tion  of  the  bulk  fluid.  That  is,  a  concentration  gradient  is 
formed.  Several  processes  like  slow  diffusion  of  the  reac¬ 
tants  in  the  electrode  pores,  diffusion  of  the  reactants  and 
products  to/from  the  electrochemical  reaction  site,  etc.  may 
contribute  to  concentration  polarization.  Concentration 
polarization  is  given  as  a  function  of  limiting  current, 
/L(v,  y)  by  the  following  expression  [2]: 


,  ,RT  f  I(x,y) 
Iconcfay)  „Fln('  IL(x,y) 


(28) 


2.3.3.  Membrane  properties  equations 

The  water  content  and  the  flux  of  water  in  the  membrane 
are  important  for  current  density  prediction  in  the  cell  [8,12] . 
The  equations  for  the  electrolyte  membrane  properties  are 
empirical  equations  and  are  obtained  by  experimental  data 
done  by  some  researchers  [8,14,15].  It  may  be  noted  that  all 
the  properties  referred  to  in  this  section  are  the  local  proper¬ 
ties  values  on  the  X-Y reaction  interface,  and  the  suffix  (x,  y ) 
has  been  taken  out  just  for  mere  simplicity.  The  activity  of 
water  (a)  on  the  membrane-electrode  interface  is  given  by 


ak 
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rw:k 


(29) 


where  subscript  k  (‘k’  =  ‘a’  or  ‘c’)  refers  to  either  the  anode 
or  cathode  side,  Xw^k  is  the  mole  fraction  of  water  at  ‘k’ 


interface,  P  the  pressure  and  P™\  is  the  saturation  vapor 
pressure  of  water  at  ‘k’  interface  and  is  given  by  [13]: 

logio(0  =  2-95  X  10“2rk  -  9.18  x  10“5rk2 

+  1.44  x  lO-7^ -2.18  (30) 

where  is  in  atm  and  Tk  the  temperature  at  the  ‘k’ 
interface  in  °C.  The  water  content  (2)  in  the  membrane  is 
given  as  a  function  of  the  activity  of  water  on  anode  side  by 
the  expression  [8]: 

2  =  0.043  +  17.81ak  —  39.85a{:  +  36 ak,  0  <  ak  <  1 

(31a) 

2  =  14+1.4(ak-  1),  1  <  <2k  <  3.  (31b) 

The  proton  (H+),  which  is  a  charged  particle,  attracts  the 
polar  water  molecule  and  drags  it  along  with  from  anode  to 
cathode  side.  This  is  called  as  the  electro-osmotic  effect. 
More  accurately,  it  is  the  hydronium  ion  [H+(H20)n]  ion 
which  passes  through  the  membrane.  The  empirical  equation 
describing  this  effect  is  [13] 

«d  =  0.002922  +  0.05/1  -  3.4  x  10“19  (32) 

where  nd  is  the  electro-osmotic  drag  coefficient.  The  diffu¬ 
sion  coefficient  of  water  ( Dw )  in  the  membrane  at  tempera¬ 
ture  Tc  (K)  is  given  by  the  expression  [13]: 

Av  =  D,  exp  ^2416  x  Qr  _  ,  T0  =  303  K  (33) 

where 

D,  =  10"10,  2  <  2 

Di  =  10"10(1 +2(2-2)),  2  <  2  <  3 

Dt  =  10-10(3  —  1.67(2  —  3)),  3  <  2  <  4.5 

Di  =  1.25  x  10"10,  2  >4.5. 

The  water  concentration  (CW;k)  at  membrane-electrode 
interface  is  given  by  [8]: 

CWik  =  ^L  (0.043  +  17.8ak  -  39.85a2  +  36a2),  ak  <  1 

44m,  dry 

(34a) 

Cw,k=^2 (14+1.4(ak- 1)),  ak  >  1  (34b) 

2kZm,dry 

where  pm? dry  is  the  density  of  dry  membrane  (2000  kg/m3 
[6])  and  Mm?dry  the  equivalent  weight  of  dry  membrane 
(1.1  kg/mol  [6]).  The  local  membrane  proton  conductivity 
((jm)  is  then  given  by  [13]: 

am=  ( 0.5 14  ^^yCwia- 0.326 

V  Pm, dry 

T0  =  303  K.  (35) 
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The  net  water  transfer  coefficient  (a)  in  the  membrane  is 
given  by  the  expression  [14]: 


n  kp  A  P  w 

t-w,m  ~ 

/% 20  _ 


(36a) 


where  Cw?m  is  the  concentration  of  water  in  the  membrane, 
kv  the  permeability  of  water  in  membrane,  /%20  the  viscosity 
of  water.  The  effect  of  third  term  (convection  term)  in 
Eq.  (36a)  is  very  small  [13],  and  is  neglected  in  this  study. 
The  modified  expression  for  net  water  transfer  coefficient  (a) 
in  the  membrane  is  then: 


codes.  All  the  equations  were  solved  and  the  parameters 
were  updated  at  the  start  of  every  iteration  using  the  adjust 
function.  UDFs  codes  were  written  in  the  C  language  and 
were  passed  as  an  input  to  the  Fluent®  software  during  the 
post-processing  step.  Because  of  the  memory  and  time 
requirements  for  solving  these  equations,  the  Fluent®  and 
the  Gambit®  software’s  were  run  on  Dell®  Precision® 
workstation  (Intel®  Xeon®  2.2  GHz,  1.5  GB  SDRAM). 


3.  Results  and  discussion 


D « 


Cw.c  Gu, 


(36b) 


2.4.  Solution  strategy 

The  model  equations  were  solved  using  the  commercial 
computational  fluid  dynamics  software  Fluent®  5.5  using 
Gambit®  2.0  as  a  preprocessor  [9].  The  flow  domain  con¬ 
sists  of  a  membrane  electrode  assemble  (ME A),  and  the 
cathode  and  anode  flow-fields  with  metal  foams.  A  hexahe¬ 
dron  volume  element  with  20  nodes  per  element,  and  map 
scheme  was  chosen  for  meshing  the  domain.  Total  number 
of  the  volume  elements  with  the  chosen  mesh  size  was 
around  0.4  million.  The  formulated  transport  equations,  viz. 
continuity,  momentum  balance,  and  species  flux  equations 
were  solved  using  the  Simple  algorithm  [9]  for  each  of  the 
control  volume  element  in  the  domain.  Thermodynamic 
equations  for  the  electrochemical  reactions  and  equations 
for  membrane  properties  were  incorporated  in  the  source 
terms  with  the  introduction  of  user  defined  functions  (UDFs) 


Simulations  were  done  for  different  permeability  levels  of 
the  metal  foam  flow-field.  The  permeability  of  the  metal 
foam  was  varied  from  10-6  to  10-12  m2.  Simulations  were 
performed  at  temperature  and  pressure  of  350  K  and 
202  kPa.  Externally  humidified  reactant  gases  were  used. 
The  anode  reactant  gas  consists  of  40  mass%  of  H2  and 
60  mass%  of  H20  vapor.  On  the  cathode  side,  the  reactant 
gas  was  humidified  air  containing  21  mass%  02,  70.5 
mass%  N2  and  8.5  mass%  H20  vapor.  The  reactants  velo¬ 
cities  at  the  anode  and  cathode  inlets  were  fixed  at  2.2  and 
4.5  m/s,  respectively.  The  cell  voltage  was  fixed  at  a  value  of 
0.6  V. 

Fig.  2  shows  the  variation  of  average  current  density  with 
the  permeability  of  the  metal  foam  used  in  the  flow-field  of 
the  bipolar/end  plate.  As  was  expected,  the  performance  of 
the  cell  increased  steadily  when  the  permeability  was 
decreased.  At  very  low  permeability  values,  comparable 
to  that  of  the  electrodes,  the  increase  in  performance  starts 
to  flatten  up.  So  it  is  clear  that  lower  permeability  values  of 
the  flow-field  will  enhance  the  fuel  cell  performance.  As 


Permeability  of 
Electrode  Material 


Permeability  =  4.4  x  10‘8  m2 


Fig.  2.  Variation  of  average  current  density  with  the  permeability  levels  of  the  metal  foam  in  the  flow-field  of  the  bipolar/end  plates. 
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Fig.  3.  Picture  of  SS-316  bipolar  plate  with  multi-parallel  flow-field 
design,  designed  and  developed  at  The  University  of  Alabama  [4]. 


already  stated,  lower  values  of  permeability  of  the  flow-field 
are  difficult  to  achieve  in  the  machined  channeled  flow-field. 
Consequently,  the  use  of  metal  foam  will  be  preferred.  We 
did  the  simulations  for  multi-parallel  design  with  channel 
in  the  flow-field  (Fig.  3),  to  compare  its  performance  with 
that  of  the  porous  metal  foam  flow-field.  The  equivalent 


permeability  of  this  channel  design  was  calculated  using  the 
expression  [5]: 


Nc{Az/2)2Axc 

L 


X 


where  Nc  is  the  number  of  channels,  Avc  the  channel  width, 
Az  the  thickness  of  the  channels  and  L  the  length  of  the 
channel.  The  value  of  the  equivalent  permeability  for  the 
multi-parallel  flow-field  channel  design  (Nc  =  24, 
Axc  =  0.0625  in.,  A z  =  0.0625  in.,  L  =  3  in.)  was  4.4  x 
10-8  m2.  The  average  current  density  obtained  for  this  case 
was  7019  A/m2  (square  mark  in  Fig.  2).  The  value  for 
average  current  density  using  the  metal  foam,  but  the  same 
permeability  value  was  6794  A/m2.  It  may  be  noticed  that 
the  average  current  density  for  the  metal  foam  is  slightly 
lower  that  that  for  the  channel  design  with  same  perme¬ 
ability.  However,  the  permeability  of  the  metal  foam  can  be 
decreased  to  improve  the  fuel  cell  performance.  But  on  the 
other  hand,  it  is  practically  difficult  to  reduce  the  perme¬ 
ability  of  the  channel  flow-field  below  a  certain  level 
(typically  10-8  m2),  due  to  difficulty  in  machining  thinner 
channels.  Also  it  is  important  to  note  that  with  use  of  metal 
foam,  a  more  uniform  distribution  of  current  density  is 
obtained.  Figs.  4  and  5  shows  that  distribution  of  local 
current  density  on  the  electrode-membrane  interface  for 
the  channel  and  metal  foam  type  flow-field.  It  is  clearly  seen 


Fig.  4.  Distribution  of  local  current  density  on  the  electrode-membrane  interface  for  the  multi-parallel  flow-field  channel  design.  Average  current 
density  =  7019  A/m2. 
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Fig.  5.  Distribution  of  local  current  density  on  the  electrode-membrane  interface  with  metal  foam  in  the  flow-field  of  the  bipolar/end  plates.  Average  current 
density  =  6794  A/m2. 


that  in  case  of  using  metal  foam  a  more  uniform  distribution 
of  local  current  density  is  observed. 

4.  Conclusions 

A  unified,  three-dimensional,  steady-state  numerical 
mass-transfer  single  cell  model  for  polymer  electrolyte 
membrane  fuel  cell  (PEMFC)  with  metal  foam  in  the 
flow-field  was  developed.  Simulations  were  performed  using 
various  permeability  levels  of  the  metal  foam  ranging  from 
10-6  to  10“ 12  m2.  It  was  observed  that  the  average  current 
density  increases  from  5943  to  8425  A/m2  when  the  perme¬ 
ability  of  the  metal  foam  flow-field  was  decreased  from  10-6 
to  10-12  m2.  However,  the  rate  of  increase  of  current  density 
decreases  when  the  permeability  values  were  decreased  to 
the  order  of  that  of  the  electrode  permeability  (  10-12  m2). 
Further,  with  the  use  of  metal  foam,  a  more  uniform  dis¬ 
tribution  of  local  current  density  was  obtained  compared  to 
that  for  the  channel  flow-field  design.  The  use  of  metal  foam 
in  the  flow-field  of  the  bipolar  plate  with  proper  permeability 
values  will  surely  enhance  the  performance  of  the  fuel  cell. 
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